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1. Introduction 

We briefly review recent developments in the mini black hole production and evaporation 
mainly based on the series of works done by Ida, Oda and Park [pi pL BL ffl]. In the low energy 
gravity scenarios such as ADD and RS-I, the CERN Large Hadronic collider (LHC) will become 
a black hole factory [|], g]. Above the TeV Planck scale, the classical production cross section of 
the (4 + n) -dimensional black hole grows geometrically o ~ s l ^ n+1 \ with a/? being the center of 
mass energy of the parton scattering. 

Once produced, black hole loses its energy or mass primarily via Hawking (thermal) radia- 
tion. The Hawking radiation goes mainly into the standard model quarks and leptons (spinors) and 
gauge bosons (vectors) localized on the brane, except for a few gravitons and Higgs boson(s). The 
quanta of Hawking radiation will have characteristic energy spectrum determined by the Hawking 
temperature and the greybody factor. The process of Hawking radiation in four dimensional rotat- 
ing black hole has been treated in detail by Teukolsky, Press, Page and others in 1970s'. In higher 
dimensions, however, it is shown that the process has quite different features. 

• Hawking temperaure 7# oc (Mfo/ ! )/M*) 1 /(" +1 ) of a (4 + n) dimensional black hole is much 
higher than 4 dimensional one with the small fundamental scale M* ~ TeV <C Mpi anc k. With 
this high temperature, the number of available degrees of freedom for Hawking radiation 
are much bigger in (4 + n) dimensions with all the standard model particles localized on the 
brane 

• The near horizon geometry of (4 + n) dimensional black hole is quite complicated. Its ge- 
ometry is different from that of a four dimensional Kerr black hole. With the highly mod- 
ified geometry in the vicinity of the event horizon, frequency dependent correction factor 
of Hawking radiation, i.e., greybody factor, is also largely modified (also see the references 



, 10, 11] as independent studies on the same topic). 



To understand the physics of those black holes, we have to understand the greybody factor of 



higher dimensional, rotating black hole [12, 13]. 



2. Generalized Teukolsky equation and greybody factor 

The induced metric on the three-brane in the (4 + n) -dimensional Myers-Perry solution [{7|] 
with a single nonzero angular momentum is given by 

A-a 2 sin 2 #, 7 lair 2 + a 2 - A) sin 2 # , , 
g = - dt z -\ - — dtdcp 

(r 2 + a 2 ) 2 -Aa 2 sin 2 # . , , „ , , 

-i— - — sin 2 tidy 2 --dr 2 -Ld$ 2 , (2.1) 



A 

where 



I = r 2 + a 2 cos 2 #, A = r 2 +a 2 - \lr x ~ n . (2.2) 
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The parameters }X and a are equivalent to the total mass M and the angular momentum J 

M= (2 + n)A 2+n ^ J = A^a ^ 
\67tG4 +n SnG4 +n 

evaluated at the spatial infinity of the (4 + n) -dimensional space-time, respectively, where A 2 + n = 
2j E (3+»)/ 2 /T((3 + n)/2) is the area of a unit (2 + «)-sphere and G\+ n is the (4 + n) -dimensional 
Newton constant of gravitation. 

1. Subtracting outgoing wave contamination at NH and separating ingoing and outgoing wave 
at FF are described. 2.Here we answer the question :what fraction of energy would be radiated into 
Hawking radiation in spin-down phase. 

2.1 Asymptotic solutions in Kerr-Newman frame 

We are given a linear, second-order equation, say 

d 2 R dR 

— + „- + I tf = 0, (24) 



where tj and t are determined in Kerr-Newman frame as: 

(s-l)A' + 2iK 

71 = A ' 

2i(Or(2s - \) - X 



n-l 



(2.5) 
(2.6) 



A = r 2 + a 2 -(r 2 H + a 2 ){^ r ) 

K = (r 2 +a 2 )w-ma. (2.7) 

The asymptotic solutions are given at NH and FF limits: 

R Nli ~ W in + W out e 2ikr * A s , (2.8) 

R W ~ Y^r 2 *- 1 + Y oxA e 2ikr * jr. (2.9) 

2.2 BC: Subtracting outgoing contamination at NH 

The solutions near the horizon r — > r# are 

Rf* = l +ai (r-r H ) + ^(r-r H ) 2 + --, 

Rout = e mr * (r ~ r H ) s (1 + h (r - r H ) + ••■) , (2.10) 
where the coefficients a,-'s and bfs are straightforward to compute: 

ax = - — , (2.11) 
1—1 

(TJo + T-Qai+Tp 

a 2 — , (2.12) 

(2.13) 



3 



Black hole 's Life at colliders 



Seong Chan Park 



where rjj and Zj are 7-th order coefficients of Taylor expansion of Tj and X, respectively. 
For 5 = 1/2 and 1, 

2/6)4,1 - A 

T_i = , (2.14) 



A 



To = A^ + 5,. 1 ^(A 1 -A 2 ), (2.15) 
Af Af 

r\-i = ^,1/2-2^;, (2-16) 



where 



T] = -^(^lAi-JToAa), (2.17) 



Ai = 2 + («-l)(l+a 2 ), (2.18) 

A 2 = l-n(?i-l)(l+a 2 )/2, (2.19) 

K = (l+a 2 )(0-am, (2.20) 

^1 = 2(0. (2.21) 



The problem is to integrate Eq-^J from purely ingoing initial conditions at r = ro out to r — > 00. 
Choosing the positive choice for s makes Y[ n stable and easily determined by an outward integration. 
However, for such an integration is unstable against contaminating the purely ingoing solution. 
We can avoid the difficulty in mid-integration, relying on a mathematical transformation of the 
equation to stabilize the solutions in the two asymptotic regimes r — > ru and r — > 00 as follows. To 
counteract the above contamination, let 

R = R-(\+ ai (r-r H )). (2.22) 

Then / satisfies the equation 

&R = g, (2.23) 

where J? = d 2 /dr 2 + r\d / dr+z wdg = —££ (1 + ai(r— r H )) = -r\a\ — T (1 +a\{r — r H )). Equa- 
tion 2.23 is now stably integrated through both asymptotic regimes, i.e., from the near horizon to 
far field regimes. 

Near the horizon, R becomes 



R(r^r H ) = -j-{r-r H f + ---. (2.24) 



2.3 Separating solution at FF 

For 5 = 1/2, 



R 



FF 



r? 
K, 



FF 



f f f 



-if 



if 



(2.25) 
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where 



cf = -t±, (2.26, 

(2.27) 



Then, R^^ becomes 



Ri/i{r - °°) - (*in + 1 - fli) +a x r + + Y out . (2.28) 

r r 

Using this expression, we can easily read out Y's without numerical difficulties. 
Finally, the greybody factor could be written as 



r-V2 = i-^w- (2-29) 



For 5=1, 



■y f* ]'* 



RZ~e 2ik ^(l + d -L + d 4 + d 4 + ■-), (2-30) 



where 



r r r- r J 



c? = -/A, (2 .31) 

_ A 2 -4fl0)(a0)-m) 
C2 " 8^ ' (2 - 32) 



Then, becomes 

y. c v e 2icor t 

R^ r _> oo) ~ (y in - fll )r+ (y ta c? - 1 +ai) + + ^out • (2.33) 

r r 

Using this expression, we can easily read out Y's without numerical difficulties. 
Finally, the greybody factor could be written as 

r^l-^M". (2.34) 



3. Hawking radiations of mini-black hole 

Schematically black hole evolution follows five successive steps as is depicted in Fig.l: the 
production phase, the balding phase, the spin-down phase, the Schwarzschild phase and the Planck 
phase. When a black hole is produced in high energy collision (production phase), the geometry 
is highly irregular, and could even be topologically non-trivial. By emitting (bulk) gravitons and 
other particles, the black hole will be settled down to a rotating black hole which is supposed to be 
well described by Myers-Perry solution in (4 + n) dimensional spacetime (balding phase). 
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The decay in spin-down and Schwarzschild phases are calculable in terms of Hawking radi- 
ation. We are interested in those phases (spin-down phase and Schwarzschild phase) and would 
answer what fraction of energy will be lost in each of these phases. 

The rate of energy (and angular momentum) loss by Hawking radiation is given as follows: 



dco(N s j. m ) , (3.1) 









where g s is the number of "massless" degrees of freedom at temperature T, namely, the number of 
degrees of freedom whose masses are smaller than T, with spin s. The expected number of particles 
of the species of spin s emitted in the mode with spheroidal harmonics /, axial angular momentum 
m is 

™»)- e ( fl ,-™n)/r_(_ 1 )2,- (3 - 2) 




Figure 1: Hawking radiation from D = lOblackhole. s= 1, a = 0.01,0.3,0.6,0.9, 1.2, 1.5. 
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4. Time evolution 



From the ratio of energy and angular momentum in eq.|3.1|, we can define a scale invariant 
function y{a s = a/r s ) as follows: 

i . N d\na, 

' r W = 5bS (41) 



2 \adM n+l 

Now we calculate the ratio of final(M^) and initial(M,) energy of black hole by integrating the 



eq.|4. l| with a s (mi) for initial angular momentum. 



M f ( /^( nnal ) v(a s ) , 

-f =Exp / da s ^- . (4.3) 



Mi \Ja s (mi) ' Cl s 

The amount of energy which is radiated in spin-down phase (0 » a s (final) ^ a s ^ a s (im)) is (M,- — 
Mf) and Mf will be also radiated in Schwarzschild phase where the angular momentum of black 
hole is zero. 

Next, let us consider the evolution of the black hole. Since time scales as r£+ 3 in (4 + n) 
dimensions 1 , it is convenient to introduce scale invariant rates for energy and angular momentum 
as follows. 

„(„) s (4.4) 

P(a.) = -C 3 ^f, (4.5) 

with these new variables y(a s ) can be written as y~ l (a s ) = f$/a(a s ) — (n + 2)/(n + 1). We also 
introduce dimensionless variables y and z to take angular momentum and mass of the hole: 

y = -lna s , (4.6) 
«=-h.£, (4.7, 

then finally we get the time variation of energy and angular momentum in terms of scale-invariant 
time parameter (t = r^" _3 (ini)f) with initial mass of the hole: 

dz a 



dy /3-a(*±f' 
dy , „ f n + 2\. n+3 — 

srw—U+i )em - (4 - 8) 



After finding the solutions z(y) and t(v) of the coupled differential equations 4.8, one can get 
y(z) and z(t), hence a s and M/Mu as a function of time. From these, one can find how other 
quantities evolve, such as the area. 



1 We can easily understand this by simply looking at the formula —dM/dt ~ AJ 4 where the surface area of horizon 
A ~ r% for brane fields and the temperature of the hole T ~ l/r s and M ~ r" +1 . 
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Figure 2: Evolution of bh in D = 10. 

Up to now we have used a unit where the size of event horizon is fixed as r/, = 1 and angular 
momentum of the hole is parameterized by (ah = a/r^). For conversion of unit, the following 
expressions are useful with a s = ah /{I + a\) x ^ n+l \ 

a(a s ) = -C\l+al)^rl^, (4.9) 

^(0.) = -«f fl (l + ^)^~, (4-10) 

where 

ln = r s M-^ = (-———) , (4.11) 

\(n + 2)a n+2 J 

,71 + 2, 1 

K n = ln{^-)~ 1 - (4.12) 

In Fig.||, black hole evolution in units of the initial mass as a function of rotation parameter a s 
for scalar(s), fermion(f), vector(v), and sum of all the standard model particles(SM) in D = 5(left) 
and D = 10 (right) are shown. The initial angular momentum parameter is fixed by a s = 0.83 
and 2.67 in D = 5 and D = 10 that are the maximal rotations allowed by the initial collision, 
respectively. The mass of the hole goes to zero before the rotation parameter goes to zero when 
only scalar emission is available. However, when all the standard model fields are turned on, the 
evolution is essentially determined by the spinor and vector radiation. It is found that a black hole 
spins down quickly at the first stage with large rotation parameter a s and the decrease of rotation 
parameter slows down as angular momentum of the hole is reduced. 

When all the standard model fields are turned on (SM), the evolution is essentially determined 
by the spinor and vector radiation. The figures show that a black hole spins down quickly at the first 
stage with large rotation parameter and the decrease of rotation parameter slows down as angular 
momentum of the hole is reduced. 

5. Summary and Discussion 

The complete description of Hawking radiation to the brane localized SM fields and the con- 
sequent time evolution of mini black hole in the context of low energy gravity scenario has been 
made. 
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We have developed analytic and numerical methods to solve the radial Teukolsky equation 
which has been generalized to the higher dimension (D = 4+n). Two main points in our numerical 
methods are as follows. First, we have imposed the proper purely-ingoing boundary condition near 
the horizon without the growing contamination of the out-going wave by extracting lower order 
terms explicitly. Second, we have developed the method to fit the in-going and out-going part from 
the numerically integrated wave solution at far field region by explicitly obtaining the next-to-next 
order expansion (or next-to-next-to-next order in vector case) of the solution. With these progress 
in numerical treatment, we can safely integrate the generalized Teukolsky equation up to very large 
r without out-going wave contamination. 

Then we have calculated all the possible modes to completely determine the radiation rate of 
the mass and angular momentum of the hole. Totally 3407 are computed explicitly, other than the 
modes which are confirmed to be negligible. A black hole tends to lose its angular momentum at 
the early stage of evolution. However the black hole still have a sizable rotating parameter after 
radiating half of its mass. More than 70% or 80% of black hole's mass is lost during the spin down 
phase. 

Now that we have completely determined the radiation and evolution of the spin-down and 
Schwarzschild phases, only remaining hurdle is the evaluation of the balding phase, which is still 
being disputed due to its non-purturbative nature, to extract the quantum gravitational information 
at the Planck phase from the experimental data at LHC. 
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